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基于 曲率 信息 的 人 工蜂 群 点 云 配 准 算法 
付 鲁 ， 陈 雷 "” 


(1. 天 津 大 学 电气 自动 化 与 信息 工程 学 院 ， 天 津 300072; 2. 天 津 商业 大 学 信息 工程 学 院 ， 天 津 300134) 


摘 要 : 针对 仅 使 用 群 智能 优化 算法 及 点 云 空间 信息 进行 点 云 配 准时 ， 优 化 过 程 寻找 两 片 点 云 对 应 点 耗 时 较 长 ， 收 
伊 速 度 较 慢 的 缺点 ， 提 出 一 种 基于 曲率 信息 的 人 工蜂 群 点 云 配 准 算法 。 算 法 根据 曲率 信息 提取 特征 点 ， 通 过 改进 人 
工蜂 群 算 法 优化 目标 函数 得 到 可 以 使 两 片 点 云 重 合 的 最 佳 变换 天 阵 。 在 种 群 优化 过 程 中 根据 曲率 信息 约束 对 应 点 寻 
找 范 围 ， 缩 小 参与 计算 点 云 的 规模 。 对 比 实验 表明 ， 与 仅 采用 随机 选 点 方法 和 使 用 点 云 空间 坐标 信息 的 配 准 算法 等 
相 比 ， 所 提出 算法 可 以 在 不 降低 配 准 精度 的 同时 ， 有 效 加 快 配 准 收敛 速度 ， 显 著 缩短 点 云 配 准 所 用 时 间 。 
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Point cloud registration based on artificial bee colony algorithm and curvature information 


Fu Kunl, Chen Lei!? 
(1. School of Electrical & Information Engineering, Tianjin University, Tianjin 300072, China; 2. School of Information 
Engineering, Tianjin University of Commerce, Tianjin 300134, China) 


Abstract: Some methods for points cloud registration only use the swarm intelligence optimization algorithm and the points’ 
spatial information. Those methods have such shortcomings as slow convergence speed and long time consuming. To 
overcome these weakness, this paper proposed a registration algorithm based on the curvature information of points cloud. 
The algorithm extracted feature points according to the curvature information. And it obtained the best transformation 
matrix to make the two point cloud coincide by the improved artificial bee colony algorithm. In the process of population 
optimization, the corresponding points were searched according to the curvature information, and the scale of point cloud 
was reduced. The experimental results show that compared with the registration algorithm only using random point selection 
method and point cloud spatial coordinate information, the proposed algorithm can effectively accelerate the convergence 
Speed and significantly shorten the registration time without reducing the registration accuracy. 
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0 ”引言 泛 应 用 于 解决 各 种 工程 优化 问题 ,文献 [4~7] 将 粒子 群 优 化 算 
本 法 用 于 工业 系统 永 磁 同步 电机 的 参数 估计 中 ， 有 效 提 高 了 系 
三 维 成 像 是 针对 某 一 实际 物体 ， 用 三 维 扫描 设备 从 不 同 。 统 性 能 。 文 献 [8] 在 语音 识别 算法 中 引入 了 粒子 群 优化 算法 ， 
角度 采集 该 物体 三 维 点 云 数据 ， 并 形成 被 测 物 体 的 三 维 完整 ” 提高 了 语音 识别 算法 的 成 功率 。 文献 [9] 在 数字 图 像 处 理 算法 
貌 ,在 三 维 成 像 领域 , 三维 点 云 的 配 准 技术 是 至 关 重 要 的 ， 中 使 用 了 蜂 群 优化 算法 ， 解 决 了 数字 图 像 水 印 技术 中 存在 的 
它 的 任务 是 找到 一 个 最 优 欧式 变换 ， 使 同一 物体 从 不 同 角 度 误 报 问题 。 文 献 [10] 将 蚁 群 优化 算法 应 用 到 骨 入 式 系统 设计 


党 


Pp 


得 到 的 有 部 分 重 欠 的 点 云 变 换 到 同一 坐标 系 下 ， 最 终 恢复 出 中 ， 成 功 实现 了 软 硬 件 的 双 路 规划 。 

被 测 物体 的 完整 形 貌 。 在 三 维 成 像 领 域 ， 由 于 群 智能 优化 算法 优异 的 全 局 优化 
维 点 云 的 配 准 过 程 实质 上 是 一 个 最 优化 过 程 。 在 求解 能 力 ， 基 于 群 智能 优化 的 点 云 配 准 技术 也 开始 逐渐 成 为 研究 

该 最 优化 问题 时 ， 传 统 的 配 准 方法 容易 受到 量化 误差 、 点 云 热点 , 并 已 表现 出 优 于 传统 方法 的 明显 优势 。 例 如 Cordon 等 

噪声 和 点 云 初始 位 置 的 影响 。 比 如 经 典 的 点 云 配 准 算法 一 一 。 人 HW 提出 了 基于 细菌 算法 的 图 像 配 准 算法 ，Santamarfa [23 提 

迭代 最 近 点 (iterative closestpoints，ICP ) 算法 中 在 待 配 准点 出 了 基于 CHC 算法 的 图 像 配 准 算法 ; Agrawal 等 人 03] 提 出 了 

云 初始 位 置 不 佳 时 往往 会 陷入 局 部 极 值 而 导致 配 准 失败 。 尽 基于 进化 刚体 对 接 算法 的 图 像 配 准 算法 等 。 


管 已 经 有 一 些 学 者 提出 了 改进 的 ICP 算法 。 如 ，Chen 等 人 口 虽然 这 些 基于 群 智能 优化 的 三 维 点 云 配 准 算法 比较 好 
在 ICP 算法 基础 上 ， 选 取 了 点 到 对 应 点 切 平面 的 距离 中 值 而 ”地 弥补 了 传统 基于 梯度 优化 的 配 准 算法 的 缺点 ， 有 效 提升 了 
不 是 对 应 点 之 间 的 距离 中 值 作为 配 准 的 目标 函数 进行 配 准 。 点 云 配 准 过 程 的 适应 性 和 配 准 精度 。 但 该 类 算法 需要 经 过 多 


文献 [3] 采 用 主 成 分 分 析 法 先进 行 粗 配 准 , 再 采用 一 种 基于 曲 ”次 迭代 进化 得 到 精准 的 配 准 结果 ,而 每 一 次 迭代 进化 过 程 中 ， 


率 特征 点 的 改进 ICP 算法 进行 细 配 准 。 然 而 这 些 改进 算法 依 。 需 针 对 种 群 个 体 计算 目标 函数 值 ， 计 算 量 较 大 。 因 此 ， 在 保 
然 存 在 配 准 易 陷入 局 部 收敛 的 缺点 。 证 点 云 配 准 精 度 的 同时 ， 减 少 迭 代 进 化 过 程 计算 量 ， 缩 短 配 
近 些 年 ， 群 智能 优化 算法 作为 新 兴 的 最 优化 方法 ， 被 准时 间 ， 就 成 了 众多 学 者 关注 的 重点 。 
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为 此 ,本文 提出 了 
bee colony algorithm, ABC ) [ 
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2 ”改进 的 人 工蜂 群 算法 


点 准 算法 。 首 先 ， 算 法 通 
过 随机 采样 和 根据 点 云 曲率 信息 采样 两 种 方式 提取 部 分 点 云 ; 人 工蜂 群 算法 是 Karaboga 等 人 模仿 蜜蜂 采 蜜 过程 提出 
> 后， 在 寻找 对 应 点 阶段 根据 点 云 曲 率 范围 的 不 同 ， 排 除 待 。 ”的 , 已 被 广泛 应 用 于 寻找 多 维 优化 问题 中 的 最 优 解 U5~!1。 该 
配 准点 云 中 曲率 范围 相差 较 大 的 部 分 点 云 ， 从 而 缩小 了 对 应 。” 算法 将 优化 问题 中 的 每 一 个 解 看 做 是 一 个 食物 源 。 每 一 个 食 
点 的 寻找 范围 ， 减 少 算法 计算 量 ， 最 后 ， 通 过 人 工蜂 群 算法 ，” 物 源 有 优 劣 之 分 ， 其 优 劣 由 各 种 因素 决定 ， 包 括 食物 源 的 地 
优化 求解 目标 函数 得 到 待 配 准点 云 之 间 的 空间 最 优 变换 。 对 ” 理 位 置 、 食 物 的 多 少 和 获取 的 难 易 程度 等 。 对 于 每 一 个 优化 
比 实验 可 以 有 效 证 明 ， 本 算法 利用 点 云 曲 率 信息 限制 了 种 群 问题， 都 有 一 个 目标 函数 来 评价 每 一 个 解 的 优 劣 。 
过 代 进化 过 程 中 寻找 对 应 点 的 范围 ， 在 保证 配 准 精度 的 前 提 整个 蜂 群 有 三 种 蜜蜂 , 分 别 是 采 蜜 蜂 、 观察 蜂 和 侦查 蜂 。 
下 ， 有 效 缩短 了 配 准 时 间 ， 提 高 了 配 准 成 功率 。 采 蜜 蜂 各 自 寻 找到 一 个 食物 源 ， 即 多 维 优化 问题 的 一 组 解 ， 
| 每 次 优化 迭 代 都 改变 这 组 解 的 部 分 信息 得 到 新 的 解 ， 通 过 比 
1 ”点 云 配 准 模型 较 目 标 函 数值 的 优 务 ， 保 留 两 者 中 最 优 的 一 个 。 然 后 采 蜜 蜂 
点 云 配 准 的 目的 是 将 通过 传感器 扫描 得 到 的 不 同 视角 将 这 个 食物 源 的 信息 带 回 蜂 集 ， 以 蜜蜂 间 特 有 的 方式 以 一 定 
的 点 云 合并 到 同 维 坐标 系 下 , 得 到 物体 的 完整 几何 模型 。 ”概率 将 它 所 携带 的 食物 源 信息 传达 给 观察 蜂 ， 观 察 蜂 同样 对 
在 空间 坐标 系 中 ， 物 体 模 型 的 刚性 变换 可 以 用 一 个 旋转 平移 接收 到 的 信息 作出 部 分 改变 得 到 新 的 解 并 择优 保留 。 此 后 经 
矩阵 表示 ， 本 文 称 这 个 矩阵 为 7 。 对 于 待 配 准 的 动态 点 云 ” 过 一 定 次 数 的 循环 ， 如 果 经 过 规定 优化 次 数 ， 观 察 蜂 没有 得 
P={py,i=42,…,N, 和 静态 点 云 Q={9),j=42,…,N,， 就 是 要 获得 到 更 好 的 食物 源 ， 则 放弃 所 携带 的 信息 ， 转 变 为 侦查 蜂 ， 并 
两 片 点 云 间 的 最 优 欧 式 变换 和 矩阵 7 了 ， 使 两 点 云 重 合 部 分 的 对 ”寻找 一 个 新 的 食物 源 。 


应 点 在 空间 华 标 系 


的 表达 式 为 
其 中 : 
在 理 


想 状 态 下 ，P 站 


Ph 完全 重合 。 该 变换 入 
数 , 分 别 为 沿 三 个 坐标 轴 的 平移 量 7.7.7 以 及 绕 3 个 坐标 有 
旋转 角 “47 。 在 不 考虑 模型 大 小 伸缩 的 情况 下 ， 变 换 矩 阵 7 
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E 阵 包含 6 个 待定 参 


在 基本 ABC 算法 中 ， 采 蜜蜂 和 观察 蜂 共 用 同一 个 搜索 


的 


(1) 


(2) 


G) 


(4) 


(5) 


与 其 在 2 中 的 对 应 点 9; 之 间 的 距离 应 当 为 0。 实 际 上 ， 它 们 


的 距离 可 以 表示 为 


d=|7(P)-Q,l. 


(6) 


但 是 由 于 测量 误差 以 及 噪声 等 因素 的 影响 ， 对 应 点 之 间 


方程 来 对 所 携带 信息 作出 改变 ， 搜 索 方程 如 下 : 
V ,=X,,+p,(X,, —X,,) (9) 
其 中 : 是 [之 间 的 随机 数 . 
对 于 解决 复杂 度 高 的 多 维 非 线性 优化 问题 ， 人 工蜂 群 算 
法 作为 群 智能 优化 算法 了 控制 参数 少 、 全 局 收敛 


的 一 种 ， 
性 好 、 易 于 实现 等 优势 ， 通 过 不 同 工 种 蜜蜂 之 间 的 合作 ， 能 
快速 准确 的 获得 全 局 最 优 解 ， 具 有 比 传统 遗传 算法 、 粒 子 群 
算法 等 群 智能 优化 算法 更 好 的 优化 性 能 。 
然而 基本 人 工蜂 群 算法 在 局 部 区 域内 的 开发 能 力 不 足 ， 
往往 导致 优化 过 程 难 以 快速 收敛 。 针 对 此 问题 ， 同 时 受 差分 
进化 算法 (DE) 19 的 启发 ， 高 卫 峰 提出 了 一 种 EABC 算法 
(enhancing artificial bee colony algorithm, EABC) 091]。 该 算 
法 采用 以 下 两 个 公式 分 别 作 为 采 蜜 蜂 和 观察 蜂 的 搜索 方程 : 
V, = + oa(X, — K+ 
BX — Xi,) (10) 


区 

BPR as) 
其 中 : x… 为 当前 群体 的 最 优 个 体 ; m1*2*i 是 从 集合 {12.… 
中 随机 选择 的 整数 ; j 是 一 个 随机 选择 的 整数 ，je{12.…,D} 。“ 
是 [0.4] 之 间 的 一 个 随机 数 ，4 为 rand*8， 这 里 B 是 一 个 均值 
为 <， 标 准 差 为 。 的 高 斯 分 布 数 。 采 蜜蜂 和 观察 蜂 分 别 对 应 
的 两 个 搜索 方程 由 于 既 有 当前 群体 的 最 优 个 体 x 的 引导 又 
有 XxX-x 的 扰动 ， 在 一 定 程度 上 保持 了 种 群 的 多 样 性 ， 可 以 
为 搜索 过 程 带 来 更 多 有 用 的 信息 ， 所 以 可 以 产生 一 个 更 有 前 
途 的 候选 解 以 此 来 提高 算法 的 性 能 ， 且 该 搜索 方程 在 一 定 程 


的 距离 无 法 达到 理想 值 0， 所 以 点 云 配 准 问题 就 转换 为 了 优 


化 问题 : 


结合 人 工 
数 "0) : 


其 中 : 


寻求 一 个 最 优 的 变换 外 
与 9 之 间 的 欧 


氏 距 离 最 小 ， 即 


蜂 条 


E(T)= min[|7(P)-©Q,|] 


6E 阵 7 使 得 所 有 对 应 点 对 T(p) 


F(T)= MedsElld ll 


f 优 化 算法 ， 将 对 应 点 距离 中 值 作为 目标 函 


度 上 平衡 了 算法 的 探索 和 开发 能 
3 ”基于 点 云 曲 率 信息 的 人 工蜂 群 配 准 算法 


点 云 配 准 过 程 一 般 存在 计算 量 大 、 时 间 复 杂 度 高 以 及 对 

应 点 寻找 不 准确 等 问题 。 针 对 以 上 问题 ， 近 年 来 有 许多 学 者 

在 点 云 的 粗 配 准 过程 中 使 用 了 点 云 的 特征 点 信息 ， 如 局 部 曲 

率 特征 信息 PIU、FPFH 特征 信息 PC、LSP 特征 信息 中 9 等。 点 

云 特征 信息 的 使 用 ， 在 有 效 减 小 参与 计算 点 云 规模 的 同时 ， 

为 对 应 点 的 匹配 提供 了 独特 的 特征 信息 。 例 如 点 云 的 不 同位 
率 


的 中 值 。 


MedsE 表示 两 片 待 配 准点 云 中 每 一 对 对 应 点 距离 平方 4 


是 曲率 大 小 一 般 差 别 较 大 ， 利 用 点 云 的 局 部 曲率 信息 ， 可 以 
效 减 小 对 应 点 误 匹 配 的 概率 。 
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于 点 云 的 特征 点 信 
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息 目前 还 主要 运用 在 点 云 粗 配 准 


应 进一步 考虑 将 特征 点 信息 和 群 智能 优化 算法 相 结 


云 配 准 算法 。 
算法 在 初始 


央 配 准 问 题 。 因 此 ， 本 文 提出 一 种 基于 曲率 信息 的 人 工 


点 云 中 提取 采样 点 时 ， 采 用 随机 提取 和 根 


布 在 符 丁 


一 个 点 经 过 


率 信息 分 段 提 取 的 方法 相 结合 , 从 而 减 小 样本 规模 。 


提取 出 的 点 云 既 包含 点 云 


率 信息 ， 又 可 以 保证 较 均 匀 地 分 


准 动态 点 云 整体 
和 峰 群 优化 过 程 中 不 断 得 到 变换 矩阵 7。 采 样 点 中 的 每 
7 变换 后 ， 仅 在 相近 曲率 范围 内 的 静态 点 云 中 寻 


记 围 内 。 
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行 采样 时 ， 从 1 到 n 中 随机 选 
一 个 数字 k， : 贡 果 动态 点 云 的 子 集 # 不 为 空 ， 且 在 静态 点 
则 从 忆 中 随机 选取 一 点 作为 


。 循 环 此 步 又， 直到 得 到 固定 个 数 的 曲率 特 
上 率 信 息 提 取 *-- 个 特征 点 ， 则 上 


云 提取 过 程 如 图 2 月 


以 合 笠 [大 


2&… 表示 根据 曲率 信息 提取 的 点 云 


-可 


找 对 应 点 。 然 后 再 计算 所 有 对 应 点 的 距离 中 值 平方 作为 蜂 群 


针对 点 云 ? 和 2 


优化 算法 的 目标 函数 ， 寻 找 最 优 变换 。 以 达到 减少 计算 量 ， 
提高 收敛 速度 的 效果 。 
言 息 的 特征 点 采样 
FP 的 每 一 个 点 ， 采 用 基于 MLS 表 务 
(moving least squares surfaces, MLS) P3] 的 方法 计算 得 到 各 


出 


基于 点 云 曲 率 


斯 | 率 KGausian 和 平 均 率 Kean o 此 得 到 加 率 : 


(12) 


(13) 


: %、 大 是 最 小 、 最 大 主 曲率 。 根 据 Chen 等 人 的 定义 ， 


以 下 公式 表示 每 一 点 的 曲率 信息 外 


k(p)+k,(p) 


1 
0, ee 
p) Se arctan (14) 


: 点 ?代表 点 云 中 的 任 
别 代表 点 ?的 最 小 主 曲率 、 

比 时 点 云 可 以 表示 为 P={Psj=02 ML ， 0={5n}sj=12 
于 点 云 模型 的 不 同 ， 


k(p)-k(p) 
E 意 一 点 ; 45) 、&(p) 和 0<5(p)<1 分 


最 大 主 曲 率 和 该 点 的 曲率 信息 。 


可 能 有 些 模型 的 点 会 大 量 吏 集 在 


率 范围 内 。 如 果 仅 仅 依靠 曲率 信息 提取 采样 点 ， 可 能 
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造成 采样 点 云 不 能 均匀 地 分 布 在 初始 点 云 的 整体 范围 内 ， 进 


会 | 


chinaXiv 


子 集 ,t=12… 


能 导致 点 云 配 准 陷入 局 部 收敛 而 失败 。 


本 文采 用 随机 采样 和 


根据 曲率 信息 采样 相 结 合 的 方式 


提取 点 云 。 点 云 中 的 每 个 点 都 


1 己 的 曲率 信息 5, 把 $s 从 0 


平均 分 为 n 段 ， 这 样 每 个 点 的 曲率 信息 就 属于 其 中 的 一 


大 一 1 


率 信 息 所 属 范围 不 同 ， 动 态 点 云 ? 可 分 为 n 个 点 云 
.no。 点 云 子 集 有 中 仅 包 含 曲 率 信息 5 满足 式 (15) 


-8<< (15) 


n n 


可 以 分 为 n 个 点 云 子 集 4 ，w=12…n。 点 云 


了 拖 集 2 中 仅 包含 率 信 息 


5 满足 式 (16) 的 点 。 


w—1.5 w+0.5 


n 


宝 2 (16) 
n 


图 1 所 示 ， 这 是 n=3 时 ， 根 据 曲率 信息 5 划分 点 云 的 


点 云 区 域 划分 


Fig.1 Region division of point cloud 


的 方式 在 人 PH 


rr 


ee 要 在 优化 过 


’ 中 的 每 一 点 在 静态 点 云 。 中 相对 
le et en ess 
秆 间 都 用 在 了 寻找 对 应 点 上 。 


的 时 间 ， 提 出 一 种 由 点 云 曲 率 信 


范围 相差 较 小 ， 所 以 在 寻找 动态 点 
， 在 静态 点 云 " 中 的 对 应 点 时 ， 无 须 在 2 的 整体 
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期 


取 
点 


征 
点 
云 


因素 导致 的 配 准 失败 ， 
提取 个 采样 点 ， 两 部 
准 的 子 点 云 。 通 过 反复 实验 
人 


分 
证 


程 


泛 忆 


云 


加 2 ” 样 点 采集 示意 图 


Fig.2 Sample collection schematic diagram 


这 种 寻找 对 应 反 


的 方法 可 以 有 效 减 少 函数 曲线 局 部 


» 寻 为 只 需 击 安 在 曲率 信 


减少 了 碍 找 对 应 点 的 四 


也 缩小 了 查找 范围 ， 从 


到 1 所 示 ， 认 为 动态 点 云 
率 信息 只 和 表态 吉 云 子 集 2 中 点 的 曲率 信息 


E 曲 率 信息 相差 较 小 的 范围 内 寻找 即 可 。 
再 计算 对 应 点 距离 平方 中 值 ， 所 得 结果 
作为 优化 算法 的 目标 函 交 


极 


息 相 


而 


子 


是 
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XIV: 


china 


录用 定稿 付 ， 蚀 ， 等 
相近 的 。 对 于 = 中 的 点 只 需 在 2 中 寻找 对 应 点 。 此 外 o 中 点 
的 曲率 信息 范围 是 包含 中 点 的 曲率 信息 范围 的 ， 这 样 就 避 
免 了 在 用 曲率 信息 分 段 边缘 的 点 不 能 准确 找到 对 应 点 。 曲 率 
信息 为 的 点 ， 其 所 属 点 云 子 集 的 序号 k 可 以 用 式 (17) 所 
未 。 

=|s/D +1 (17) 


其 中 : L 4 表示 向 下 取 整 。 为 了 加 快 搜索 速度 


度 ， 可 以 对 静态 


点 云 的 每 一 个 子 集 o% 建立 kd-tree 数据 
以 下 列 出 了 优化 算法 目标 函数 ”2) 的 求解 流程 : 
for i=1,2,…,N, 
a) 使 用 优 
b) 根 据 式 17 计算 得 到 >” 的 k 值 ; 

c) 在 %, 中 寻找 bp 的 对 应 点 4 ; 
d) 计 算 7Cz) 与 4 的 欧式 距离 4 =|76p)-gl; 
End 


和 得 到 


标 函 数值 。 


e) 由 公式 wC) = MedsElz 上 计 


4 ”实验 结果 及 分 析 
4.1 三 维 点 云 数 据 
配 准 实验 所 用 三 维 点 


Tele、Bird、Angel 以 及 frog 点 云 模 型 。 四 种 


寺 构 进行 索引 。 


化 过 程 得 到 的 变换 矩阵 7 对 pb, 进行 旋转 平移 ,得 到 点 7(p) ; 


i 云 模型 分 别 为 SAMPL 点 云 库 中 的 


点 云 模 型 都 选取 


0 度 视角 和 40 度 视角 的 两 片 点 云 。 各 模型 0 
个 数 分 别 为 6841、9 115、 
的 数量 分 别 为 5612、5 872、 
云 静止 不 动 ， 另 一 片 点 云 通过 随机 旋转 平移 
点 云 进 行 配 ; 
4.2 配 准 实验 设计 


度 视角 所 含 点 的 


14 089 和 9 640，40 度 视 角 所 含 点 
13 190 和 7 650。 配 准时 一 片 点 


变换 后 作为 动态 


为 了 分 析 比 较 本 文 算法 的 配 准 性 能 ， 进 
实验 1: 把 使 用 了 点 云 曲率 信息 的 改进 


行 两 组 实验 。 
人 工蜂 群 配 准 算 


法 (以 下 简称 本 文 算法 ) 和 仅 使 用 点 云 空间 多 
人 工蜂 群 配 准 算法 以 下 简称 空间 算法 ) 作 
为 实验 1 中 两 种 算法 所 用 策略 对 比 。 

实验 2: 把 本 文 算法 和 其 他 两 种 较为 新 
优化 算法 的 点 云 配 准 方法 进行 对 比 。 


坐 示 信 息 En 的 改进 
对 比 。 表 1 所 示 


颖 的 基于 群 智 能 


表 1 实验 1 算法 策略 对 比 
Table 1 Algorithm strategy comparison of experiment 1 

算法 优化 算法 样 点 采集 标 函 数 

率 特征 点 + 率 约束 + 

Re 随机 采样 距离 中 值 

空间 算法 EABC 随机 采样 距离 中 值 

4.3 参数 设置 

实验 用 计算 机 硬件 配置 为 Intel Core (TM) i7-4790 处 


YH 


里 器 ， 主 频 3.6 GHz， 内 存 16 GB 。 操 作 系统 为 Windows 8， 


所 有 算法 均 在 matlab R2014a 软件 下 编程 实现 。 
本 文 使 用 的 蜂 群 优化 算法 种 群 数 目 过 多 会 降低 点 云 配 
准 效率 ， 而 种 群 数目 过 少 则 会 导致 配 准 可 能 无 法 成 功 。 由 预 


实验 得 到 一 组 本 文 算法 的 控制 参数 。 


为 宜 , 其 中 采 蜜蜂 和 观察 蜂 各 占 一 半 。 最 大 党 
高 斯 分 布 均值 “=-03 ,标准 差 5-03 ，4=1。 配 } 


模型 0 度 视角 的 点 云 作为 动态 点 云 ”， 从 ? 


EABC 中 种 群 数 量 取 20 


试 次 数 limit=50， 
过 程 中 选取 各 
中 选取 300 个 点 


: 基于 曲率 信息 的 人 工蜂 群 点 云 配 准 算法 
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表 2 EABC 算法 参数 选取 
Table 2 Parameter selection of EABC algorithm 
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参数 类 别 参数 值 
种 群 数量 20 
limit 50 
采 蜜 蜂 (观察 蜂 ) 数量 10 
旋转 角度 范围 [0,360] 
平移 范围 /mm [-40mm,40mm] 
高 斯 分 布 均值 * 0.3 
标准 差 5 0.3 
A 1 
4.4 实验 1 结果 分 析 
经 过 多 次 实验 对 比 得 出 ， 实 验 1 中 的 两 种 算法 大 约 在 


160s 左右 都 可 以 收敛， 
相同 的 时 间 内 ， 每 隔 
用 向 量 ? 对 完整 动态 点 云 中 的 点 
换 ， 得 到 点 云 ro) 。 


因此 约束 实验 运行 时 间 为 200 s。 
10s 记录 下 蜂 群 中 的 最 优 变换 向 量 ”。 


在 


，i=122w 做 旋转 平移 变 


计算 Tw 与 静态 点 云 《 中 对 
r=1,2,.…,20 。 其 中 ， 
的 值 ， 共 记录 20 次 。 
的 距离 平方 中 值 。 最 后 用 4， 
观察 算法 收敛 对 比 效果 。 


Li 


r 代表 从 选 代 优化 3 


应 点 的 距离 平方 中 值 = ， 
二 始 每 隔 10 s 记录 一 
/表示 算法 运行 bxrs 时 ， 所 
,20 构造 算法 的 收敛 曲线 ， 
此 外 ， 针 对 各 点 云 数据 模型 ， 


次 上 
有 对 应 点 


分 别 


设 定 判定 配 准 是 否 成 功 的 目标 函数 限定 病 值 "。 


A Tt 出 当 
的 阔 值 。 


1 定 配 准 成 功 ， 


表 3 判定 各 模型 配 准 是 


Threshold value to determine whether registration of each 


Table 3 


否则 视 为 失败 。 表 3 为 针 


否 成 功 的 闪 值 


model is successful or not 


标 函 数值 
对 各 模型 


模型 


TtT(mm) 


Bird 
Tele 
Angle 
Frog 


0.60 
0.40 
0.70 
0.60 


姑 为 使 用 群 智能 优化 算法 进行 点 


\ 云 配 准 


在 较 小 概率 会 发 生 算法 局 部 收敛 的 
据 对 平均 数据 
配 准 失败 时 的 数据 ， 仅 取 配 准 


存在 随机 性 ， 存 


情况 ， 失 败 的 点 云 配 准 妆 
时 影响 较 大 ,所 以 两 种 算法 各 自 运行 
成 功 时 的 平均 数据 对 算法 进行 


由 于 点 云 初始 位 置 不 确定 
制 收敛 曲线 的 < 的 值 可 能 
的 收敛 
来 绘制 收 全 


， 大 


线 。 图 3 为 各 点 云 模 


较 可 以 得 出 ， 
本 文 算法 一 般 在 
敛 速度 较 慢 ， 


两 种 算法 都 可 
车 100s 之 内 可 


以 达 


引入 了 点 云 曲 率 信 息 ， 通 过 
从 而 大 大 加 快 算法 收敛 速度 。 

表 4 记录 了 两 种 对 比 算法 30 
成 功 的 次 数 中 求 取 目标 函数 的 平均 


率 的 


与 最 终 收敛 时 的 
线 不 便于 观察 收敛 效果 ， 所 以 使 用 50s 
线 。 

将 实验 所 得 各 点 云 模型 的 目标 函数 值 与 限定 阔 值 * 作 比 
以 达到 收敛 ,观察 图 3 


一 般 在 150s 左右 才 可 以 收敛 。 
本 文 算 法 收敛 速度 更 快 。 这 是 因为 在 本 文 算 法 的 目标 函 
函数 复杂 度 ， 


验 运行 开始 阶段 ， 
值 相差 较 大 ， 绘 制 


30 次 ， 


型 的 收敛 


到 收敛 ， 而 空 


限制 降 


低 了 


相同 时 间 


之 后 的 < 值 


可 以 看 出 ， 
zs 间 算 法 则 收 
内 ， 
数 中 


次 运行 中 的 成 功 次 数 ， 在 


值 以 及 200 s 


的 平均 值 。 可 以 看 出 本 文 算法 与 空 


间 算 法 相 比 ， 


行 30 次 仅仅 在 angel 模型 实验 中 


作为 采样 点 | 


ed 


] 于 配 准 过 程 。 各 模型 动态 点 云 ? 在 用 于 配 准 前 ， 
其 空间 坐标 在 [-40mm,40mm] 间 随机 平移 ， 在 o2n 的 角度 范 
7 别 绕 X,Y，Z 轴 随 机 旋转 。 表 2 为 EABC 优化 算法 3 


[TT 


而 最 终 达到 的 精度 较 之 空间 算法 并 
时 间 内 ， 本 文 算 法 迭代 次 数 明显 多 
这 说 明 本 文 算法 降低 了 每 一 次 迭代 


次 失败 ， 


内 连 代 次 数 
本 文 算法 运 


成 功率 更 高 ， 


没有 降低 。 


,本 文 算法 的 和 迭代 次 数 更 多 , 从 而 


使 算法 更 


同时 ， 在 相同 
于 空间 算法 的 迭代 次 数 。 
的 过 程 计 算 量 。 
快 地 达到 收敛 。 


相同 时 间 
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本 文 算法 | 


示 文 站 法] 
空间 算法 | -这 间 昌 法 | 


目标 函数 值 
目标 函数 值 


100 50 100 1 
运行 时 间 ( 秒 ) 运行 时 间 ( 秒 ) 


(a)Bird (b)Tele 


(0) Angel (d) Frog 
- a 


目标 表 数 信 
目标 本 数值 


100 150 100 150 
运行 时 间 ( 秒 ) 运行 时 间 (种 ) 


(c)Angel (d)Frog 

图 3 实验 1 中 各 模型 收敛 曲线 

Fig.3 Convergence curves of models in Experiment 1 
表 4 实验 1 结果 对 比 


Table 4 Comparison of results in experiment 1 


点 云 模型 配 准 算法 成 功 次 数 平均 目标 函数 值 (mm) 迭代 次 数 
. 本 文 算法 30 0.2784 6614 
Bird ok 
空间 算法 29 0.3057 3891 
是 本 文 算 法 30 0.2854 6974 
ele 
空间 算法 29 0.2698 4227 
本 文 算法 30 0.4236 6283 
Angel on 
空间 算法 26 0.4508 3653 
本 文 算法 29 0.3804 5766 
Frog i 
空间 算法 24 0.4065 3627 
对 于 四 种 模型 ， 采 用 本 文 算法 在 200s 内 迭代 次 数 分 别 


为 6614 次 、6 974 次 、6 283 次 、5 766 次 ， 采 用 空间 算法 的 
和 迭代 次 数 分 别 为 3 891 次 、4 227 次 、3 653 次 、3 627 次 。 可 
以 看 出 本 文 算法 与 空间 算法 相 比 ， 在 相同 时 间 内 ， 迭 代 次 数 
更 多 。 这 说 明 采 用 本 文 算法 可 以 有 效 减 少 每 一 次 迭代 进化 过 
程 的 计算 量 ， 因 而 相同 时 间 内 可 以 有 更 多 的 迭代 次 数 ， 进 而 
证 明 在 达到 相同 配 准 精 度 标准 时 , 本文 算法 计算 量 大 大 减少 。 
图 4 为 使 用 本 文 算 法 得 到 的 配 准 效果 。 从 图 中 可 以 看 出 
本 文 算法 配 准 效果 良好 。 


(a)Bird (b)Tele 


(c)Angel 


(d)Frog 
图 4 本文 算法 配 准 效果 
Fig.4 Registration effected map of presented algorithm 


付 锟 ， 等 : 基于 曲率 信息 的 人 工蜂 群 点 云 配 准 算法 
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4.5 实验 2 结果 分 析 

实验 2 中 ， 将 本 文 算 法 与 2016 年 提出 的 基于 改进 ISS 
特征 点 与 人 工蜂 群 算法 的 点 云 配 准 方法 (ABC2016) 中 和 
2014 年 提出 的 基于 侵略 模型 的 自 适 应 分 布 式 差分 进化 点 云 
配 准 算法 (DE2014) R51 作对 比 。 实 验 2 中 本 文 算法 采用 和 实 
验 1 中 相同 的 实验 参数 、 实 验 平 台 、 点 云 数 据 模型 、 目 标 评 
价 函 数 和 相同 的 实验 设计 。 其 余 两 种 算法 所 需 主要 优化 过 程 
参数 见 表 5。 


表 5 实验 2 中 各 算法 相关 参数 


Table. 5 Parameters of each algorithm in experiment 2 


算法 名 称 优化 过 程 参数 
ABC2016 种 群 数量 20，limit=100 
DE2014 种 群 数量 16， 种 群 规模 10， =* 


同 实 验 1, 每 种 算法 运行 30 次 , 每 隔 10 s 记录 用 当前 最 
优 解 所 对 应 的 目标 函数 值 ， 并 以 此 绘制 各 算法 的 收敛 曲线 。 
各 模型 收敛 曲线 见 图 5 


(a) Bi (5) Tale 
0 


本 文 算法 | 示 文 算法 
ABc20e | 


目标 条 数值 
目标 函数 值 


100 150 100 150 
运行 时 间 ( 秒 ) 运行 时 间 (种 ) 
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图 5 实验 2 中 各 模型 收敛 曲线 
Fig.5 Convergence curves of models in Experiment 2 
表 5 记录 三 种 算法 30 次 运行 中 的 成 功 次 数 ， 最 终 的 平 
均 目标 评价 函数 精度 。 
表 5 实验 2 结果 对 比 


Table 5 Comparison of results in experiment 2 


点 云 模 型 配 准 算法 成 功 次 数 平均 目标 函数 值 
本 文 算法 30 0.2637 
Bird ABC2016 25 0.5094 
DE2014 26 0.2936 
本 文 算法 30 0.2865 
Tele ABC2016 29 0.3952 
DE2014 28 0.2391 
本 文 算法 30 0.4359 
Angel ABC2016 28 0.5903 
DE2014 26 0.6248 
本 文 算法 29 0.3710 
Frog ABC2016 26 0.5739 
DE2014 24 0.3262 


与 实验 1 不 同 的 是 ， 实 验 2 未 统计 相同 时 间 内 各 算法 和 迭 
代 次 数 。 这 是 因为 各 配 准 算法 所 用 优化 方法 不 同 ， 种 群 粒子 
数目 不 同 ,导致 无 法 比较 同一 时 间 内 各 种 群 的 迭代 进化 次 数 。 
观察 图 5 可 知 ， 使 用 了 点 云 曲 率 信 息 的 本 文 算法 与 另外 
两 种 算法 相 比 较 ， 可 以 更 快 地 达到 收敛 。 同 时 由 表 5 可 以 看 
出 ， 与 另 两 种 算法 相 比 ， 本 文 算 法 并 未 降低 配 准 精度 。 此 夕 


其 他 两 种 算法 的 失败 次 数 要 多 于 本 文 算法 ， 证 明了 本 文 算法 
有 更 好 的 稳定 性 。 


5 ”结束 语 


本 文通 过 对 选取 采样 点 以 及 寻找 对 应 点 过 程 的 改进 ， 
出 了 一 种 基于 曲率 信息 的 点 云 配 准 方法 。 首 先 ， 在 提取 点 
预 处 理 阶 段 ， 根 据点 云 曲 率 信息 采样 与 随机 采样 相 结合 的 
式 各 提取 一 部 分 点 。 优 化 过 程 中 根据 点 云 曲率 信息 所 属 范 
的 不 同 减 小 寻找 对 应 点 的 范围 ， 从 而 减少 寻找 对 应 点 的 计 
量 。 最 后 用 改进 的 人 工蜂 群 EABC 算法 找到 点 运 之 间 的 最 
变换 向 量 。 通 过 实验 对 比 证 明 ， 本 文 算法 可 以 达到 良好 的 配 
准 效果 。 此 外 与 仅 使 用 点 云 空间 信息 的 蜂 群 优化 配 准 算法 ， 
以 及 其 他 两 种 算法 相 比 ， 本 文 算 法 可 以 在 保证 配 准 精度 的 同 
时 更 快 达到 收敛 。 同 时 在 相同 运行 次 数 中 ， 本 文 算法 配 准 成 
功率 更 高 。 与 传统 群 智能 优化 配 准 算法 相 比 ， 本 文 算 法 有 较 
强 的 竞争 力 。 
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